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We present a study of /Ujujuju, eeee, and /u/xee events using 1 fb _1 of data collected with the 
DO detector at the Fermilab Tevatron pp Collider at y/s = 1.96 TeV. Requiring the lepton pair 
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masses to be greater than 30 GeV, we observe one event, consistent with the expected background 
of 0.13 ± 0.03 events and with the predicted standard model ZZ and Z'y* production of 1.71 ± 0.15 
events. We set an upper limit on the ZZ and Z7* cross section of 4.4 pb at the 95% C.L. We also 
derive limits on anomalous neutral trilinear ZZZ and ZZ^* gauge couplings. The one-parameter 
95% C.L. coupling limits with a form factor scale A = 1.2 TeV are -0.28 < f£ < 0.28, -0.31 < 
/Jo < 0.29, -0.26 < /2o < 0.26, and -0.30 < f? M < 0.28. 

PACS numbers: 12.15.Ji, 13.40.Em, 13.85.Qk 



The standard model (SM) makes precise predictions 
for the couplings between gauge bosons based on the 
non-Abelian symmetries of the model. These predictions 
can be tested by studying di-boson production (WW, 
WZ, ZZ, Zj, and Wj) at particle colliders. ZZ produc- 
tion is predicted to have the smallest cross section among 
the di-boson processes. Because the decay channels with 
the smallest expected background also have the smallest 
branching fractions, it has not been observed at a hadron 
collider. Nevertheless, the final states with small back- 
grounds may provide the best opportunity to observe the 
effects of physics beyond the SM. In addition to the pro- 
duction of new particles that could decay into ZZ , vari- 
ous extensions of the SM predict large anomalous values 
of the trilinear couplings [lj ZZZ and ZZ^f* that would 
result in higher cross sections than the SM prediction. 
The direct ZZ^j* and ZZZ couplings are zero in the SM. 
Consequently, an observation of an enhancement of the 
cross section would indicate physics beyond the SM. 

Previous studies of ZZ production were made at the 
LEP electron-positron collider. The combined LEP re- 
sults are available in Ref. 0). All results are consistent 
with SM predictions. The CDF Collaboration searched 
for ZZ and W Z production in pp collisions at centcr- 
of-mass energy yfs — 1.96 TeV with the result that 
cr(ZZ) + a{WZ) < 15.2 pb at 95% C.L. [3] The pre- 
dicted ZZ production cross section in the SM is 1.6 ±0.1 
pb at the Tevatron Collider energy [l], H| • 

In this Letter, we present a search for ZZ and Zj* pro- 
duction and a search for anomalous trilinear ZZj* and 
ZZZ couplings at the Fermilab Tevatron Collider with 
the DO detector. We follow the framework of Ref. [lj], 
where general ZZV (V = Z, 7) gauge and Lorentz- 
invariant interactions are considered. Such ZZV cou- 
plings can be parameterized by two CP- violating (f¥) 
and two CP-conserving (f§ ) complex parameters. Ad- 
ditional anomalous couplings 0] can contribute when 
the Z bosons are off-shell. However, these couplings 
are highly suppressed near the Z boson resonance. We 
note that the ZZV couplings are distinct [l| from the 
Z^fV couplings probed in Z^j production in e + e~ and 
hadronic collisions. Partial-wave unitarity is ensured by 
using a form-factor parameterization that causes the cou- 
pling to vanish at high parton center-of-mass energy \fi: 
f¥ = fV/(l + s/A 2 ) n . Here, A is a form-factor scale, f% 
are the low-energy approximations of the couplings, and 
n is the form- factor power. In accordance with Ref. [H, 



we set n — 3 for all cases. The form-factor scale A is 
selected so that limits are within the values provided by 
unitarity at Tevatron Collider energies. 

We search for ZZ and Zj* production with a final 
state signature that consists of four leptons: two pairs of 
either electrons or muons. The electron and muon pairs 
can be produced either by the decay of an on-shcll Z bo- 
son or via a virtual Z boson or photon. We studied three 
final states: four muon (/z^t/i^t), four electron (eeee), and 
two muons and two electrons (/x/nee). 

Data used in this analysis were collected with the 
DO detector at the Fermilab Tevatron pp collider at 
y/s = 1.96 TeV between October 2002 and February 
2006. The integrated luminosities [7| of the three chan- 
nels are approximately 1 fb _1 each and are shown in Ta- 
ble D 

The DO detector Q is a multi-purpose detector de- 
signed to operate at high luminosity. The main com- 
ponents of the detector are an inner tracker, a liquid- 
argon/uranium calorimeter, and a muon system. The in- 
ner tracker consists of a silicon microstrip tracker (SMT) 
and a central fiber tracker (CFT) operating in a 2 Tesla 
solenoidal magnetic field. The SMT has coverage to 
pseudorapidity [9j \-q\ sa 3.0, while the CFT has cov- 
erage to \i]\ « 1.8. The calorimeter system is outside 
the solenoid and has three cryostats, one for each of the 
two endcap calorimeters (EC) and one for the central 
calorimeter (CC). The CC covers the region |7y| < 1.1, 
while endcap calorimeters extend the coverage to ~ 4. 
The calorimeters are segmented along the shower direc- 
tion with four layers forming the electromagnetic section 
(EM) and an additional three to five layers forming the 
hadronic section. This allows for electron-pion separation 
based on longitudinal and transverse shower shape. The 
muon system is outside of the calorimeters and consists 
of 1.8 Tesla iron toroid magnets with tracking chambers 
and scintillator counters mounted both inside and outside 
the toroid iron. The muon system extends to «2. 

The DO detector utilizes a three-level (LI, L2, and L3) 
trigger system. Electrons candidates are required to have 
deposited energy in the EM section of the calorimeter at 
LI and L2, and then to satisfy selection criteria on shower 
shape and the fraction of energy in the EM calorimeter 
at L3. Muon candidates are required to have hits in the 
muon scintillation counters and a match with a track in 
the LI track system. In a fraction of the data, muons are 
also required to have hits in the chambers of the muon 
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system. At L2, muon candidates are required to have 
track segments in the muon tracking detectors, and a 
minimum transverse momentum (j>t) is required. At L3, 
some muons are also required to have a match with a 
track from the inner tracker. 

Combinations of single-electron and di-electron, as well 
as single- muon triggers are used in this analysis. These 
triggers have varying pt and quality requirements on the 
lcptons. The trigger efficiency for events with four high- 
er lcptons that satisfy all other selection requirements 
is approximately 99%. In order to increase the collec- 
tion efficiency, we do not require pairs of leptons to have 
opposite electric charge. 

For the fJ.fJ.fifJ, channel, each muon is required to satisfy 
quality criteria based on scintillator and wire chamber in- 
formation from the muon system and to have a matching 
track in the inner tracker. The track matched to the 
muon must have px > 15 GeV. Additionally, muons that 
are only identified in the muon detector layers before the 
toroid are required to have less than 2.5 GeV of total 
transverse energy deposited in the calorimeter in an an- 
nulus between R = 0.1 and R = 0.4 centered around the 
track matched to the muon [10]. Finally, the distance in 
the transverse plane of closest approach to the beamspot 
for the track matched to the muon must be less than 
0.02 cm for tracks with SMT hits and less than 0.2 cm for 
tracks without SMT hits. This reduces the background 
from muons that do not originate from the primary ver- 
tex such as those from b decays and cosmic rays. The 
maximum distance between the muon track vertices for 
all muon pairs along the beam axis is required to be less 
than 3 cm. This reduces backgrounds from pairs of cos- 
mic ray muons and from beam halo. Since the charge of 
the muons is not considered, three possible Z/7* com- 
binations can be formed for each fifififi event. Selected 
events are required to have the invariant masses of both 
muon pairs above 30 GeV for at least one of the combi- 
nations. 

For the eeee channel, each electron is required to have 
transverse energy (Et) greater than 15 GeV, and to have 
\rj\ < 1.1 or 1.5 < \rj\ < 3.2, and to be isolated from other 
energy clusters. Electrons are also required to satisfy 
identification criteria based on multivariate discrimina- 
tors derived from calorimeter shower shape and tracking 
variables. Since the calorimeter fiducial region covers re- 
gions where there is little or no tracking coverage, only 
three of the four electrons are required to have an asso- 
ciated inner track. Electrons without a track match are 
required to satisfy more stringent shower shape require- 
ments. Similar to the fifififi channel, events are required 
to have the invariant masses of both electron pairs to 
be above 30 GeV for at least one of the three possible 
combinations. 

For the fifiee channel, muons are required to satisfy the 
same single muon selection criteria as in the fJ.fJ.fifJ. chan- 
nel and electrons the same transverse energy and rj cri- 



teria as in the eeee channel. Both electrons are required 
to satisfy the multivariate discriminator based on shower 
shape and parameters of the matching inner track. In ad- 
dition, electrons and muons are required to be separated 
by R > 0.2 to reduce backgrounds from muons that have 
radiated photons spatially coincident with a track. The 
invariant masses of the muon pair and the electron pair 
are required to be greater than 30 GeV. 

We determine the total acceptance using a combina- 
tion of information from MC and the data. We deter- 
mine single electron and muon identification efficiencies 
directly from data ll| , which contains high px electrons 
and muons from more than 100,000 single Z decays in 
each channel. Efficiencies are parameterized as functions 
of the relevant variables for the tracking, calorimeter and 
muon systems such as the position of the interaction 
along the beam direction and the rj of the electron or 
muon. The acceptance for all channels is then determined 



using ZZ and Z7* events generated with pythia [12 1 
and a parameterized simulation of the detector. The mo- 
menta of the muons and electrons are fluctuated in the 
Monte Carlo (MC) using angular and energy resolutions 
determined from the data, and the measured single elec- 
tron and muon identification efficiencies are taken into 
account in calculating the acceptance. The number of 
observed events, the acceptance, and the expected signal 
calculated assuming the cross section is 1.6 pb are listed 
in Table |T] for each channel. Systematic uncertainties in 
acceptance due to momentum and energy scale calibra- 
tions, angular resolution, lepton identification variation 
with luminosity as well as other effects were included. 
Taking into account correlations between these uncer- 
tainties, a total of 1.71 ± 0.15 ZZ and Z7* events is 
expected. 

The main background sources are "Z + multi-jet" 
events, events where a top-antitop (tt) quark pair is 
produced, and "combinatoric events" , which are four- 
lcpton events that survive because mis-paired leptons 
cause them to satisfy the dilepton invariant mass selec- 
tion criteria when they would not have otherwise. 

For a "Z + multi-jet" event to be a ZZ candidate, 
the jets are either mis-identified as electrons or muons or 
contain real electrons or muons from in-flight decays of 
pions, kaons, or heavy quarks. Though Z + jets (and 7* 
+ jets) events are the primary source of this background, 
it also includes those events where the Z boson or 7* is 
reconstructed from one or more "fake" leptons. The "Z 
+ multi-jet" event background is measured by first de- 
termining the probability for a jet in data to produce an 
electron or muon that satisfies the identification criteria. 
The probability for a jet to mimic a muon was measured 
in two-jet events selected via jet triggers. Muons that 
satisfy the selection criteria, have pr > 15 GeV, and are 
found within R < 0.2 of the lower-£r jet are considered 
"fake" . The probability for a jet to mimic a muon is pa- 
rameterized in terms of the tag-jet Et- It is 10 -4 for jets 
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Decay Channel 


Luminosity (pb L ) 


Acceptance 


Expected Signal 


Observed Events 


fifififi 


944 ± 58 


0.27 ±0.02 


0.46 ± 0.05 





eeee 


1070 ± 65 


0.23 ±0.01 


0.44 ±0.03 





fifiee 


1020 ± 62 


0.22 ±0.02 


0.81 ±0.09 


1 



TABLE I: The integrated luminosity, acceptance (including lepton identification efficiencies), expected number of signal can- 
didates, and the number of observed events in the three decay channels for the ZZ and Z'y* cross section analysis. A total of 
1.71 ±0.15 ZZ and Zy* events is expected. 



with Et ~ 15 GeV and increases approximately linearly 
with E T to 1(T 2 for jets with E T » 150 GeV. The prob- 
ability for a jet to mimic an electron was measured using 
a procedure similar to that described for muons. It is 
parameterized in terms of r\ and Et to account for dif- 
ferences in the EC and CC and is typically 10 -4 (10~ 2 ) 
per jet for those electrons with (without) a track match. 
A systematic uncertainty of 30% is assigned to account 
for variations in the probabilities as a result of changes 
of the lepton identification criteria that were performed 
as cross checks. The probabilities for jets to be misiden- 
tified as muons are then applied to the jets in our fifi + 
jets data to determine the background to fifififi. Simi- 
larly appropriate probabilities are applied to lit + jet(s) 
events to determine the background to eeee and fifiee. 
There, because we started from a three lepton sample, 
we correctly account for events with two real leptons, a 
photon misidentified as an electron, and a jet misidenti- 
fied as either an electron or muon. 

The background from tt events is determined from 
PYTHIA MC using a detailed GEANT-based [13| detec- 
tor simulation program and the same reconstruction pro- 
gram that was used for the data. 

The combinatoric background occurs only in the four- 
electron and four- muon decay channels. While it could 
be reduced by ~ 1/3 by requiring leptons which form Z 
bosons to have opposite signs, that would have resulted 
in a loss of signal efficiency for the high-pr leptons that 
could result from anomalous couplings. This background 
was estimated using MC simulation. It is 0.016 ± 0.003 
(0.015 ± 0.003) events in the /i/x/i/x (eeee) channel; the 
uncertainty in the background comes from uncertainty 
in the lepton px resolution. 

Table ILT1 displays the contributions of all non-negligible 
backgrounds and a summary for each decay channel. The 
total expected number of candidates from background 
sources is 0.13 ± 0.03 events. 

One event is observed in data, consistent with the 
SM prediction of 1.71 ± 0.15 events plus background of 
0.13 ± 0.03 events. We set an upper limit of 4.4 pb at 
the 95% C.L. on the cross section for pp — > ZZ + X and 
Z7* + X, where dilepton pair masses are greater than 30 
GeV. 

Because we do not observe an excess of events, we set 
limits on anomalous trilinear couplings by comparing the 
number of observed candidates with the predicted back- 



Background 




eeee 


fifiee 


Z ± multi-jet 


0.004 ±0.001 


0.065 ±0.021 


0.007 ±0.002 


tt 


0.010 ±0.005 




0.006 ± 0.003 


Combinatoric 


0.016 ±0.003 


0.015 ±0.003 




Beam Halo 


0.003 ± 0.001 






Total 


0.033 ± 0.006 


0.080 ± 0.021 


0.013 ±0.004 



TABLE II: Contributions of non-negligible backgrounds to 
the expected number of candidates in the three decay channels 
for the cross section analysis. The total expected background 
is 0.13 ± 0.03 events. 



ground and the expected sum of ZZ events from a MC 
program [l[ using anomalous ZZ"/* and ZZZ couplings. 
We produced grids of MC samples by varying two of the 
anomalous couplings at a time. There are typically 5000 
events generated at each grid point. These events are 
processed through the same detector and reconstruction 
simulation as was used in the cross section analysis. Be- 
cause on-shell Z bosons dominate the contributions of 
the ZZV couplings in the MC and because off-shell Z 
bosons are enhanced by couplings [B[ not implemented 
in the MC, we required the dilepton invariant mass to 
be > 50 (70) GeV for dimuons (dielcctrons). This selec- 
tion criterion is set by the resolution of the dilepton mass 
measurement and it removed the single fifiee candidate 
event. 

The number of events expected for each choice of 
anomalous couplings was used to form a likelihood Q 
for that point. Poisson probabilities were used for the 
expected signal plus background. These are convo- 
luted with Gaussian uncertainties on the acceptance, 
background, and luminosity. One-dimensional (ID) and 
two-dimensional (2D) limits on anomalous couplings are 
formed by finding the coupling values with likelihood of 
1.92 and 3.00 units greater than the minimum. Limits 
are determined using A = 1.2 TeV. The 95% C.L. ID 
limits are -0.28 < fg < 0.28, -0.26 < f2 < 0.26, 
-0.31 < < 0.29, and -0.30 < / 5 7 < 0.28. The 95% 
C.L. 2D contours / 4 7 vs. fg, fj vs. / 5 7 , fg vs. /|>, 
and /gg vs. fg are shown in Figure [T] In the four 2D 
cases the other two couplings are assumed to be zero. 

In summary, we report results from a search for ZZ 
and Zj* production using the eeee, /x/i/i/i and fifiee de- 
cay signatures. We analyzed 1 fb _1 of data collected with 
the DO detector at the Fermilab Tevatron pp Collider at 
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FIG. 1: Limits on anomalous couplings for A = 1.2 TeV: (a) 
f2o vs. fi , (b) f2 vs. f7 M , (c) fi vs. f§ , and (d) f7 M 
vs. fso, assuming in each case that the other two couplings 
are zero. The inner and outer curves are the 95% C.L. two- 
degree of freedom exclusion contour and the constraint from 
the unitarity condition, respectively. The inner crosshairs are 
the 95% C.L. one-degree of freedom exclusion limits. 

\/s = 1.96 TeV. Requiring each lepton to have transverse 
momentum greater than 15 GeV, and the dilepton pair 
masses to be greater than 30 GeV, we observe one event 
with an expected SM background of 0.13 ± 0.03 events. 
The one observed event is consistent both with back- 
ground and with predicted SM ZZ and Z"f* production 
of 1.71 ± 0.15 events. We set an upper limit of 4.4 pb at 
the 95% C.L. on the cross section for pp — > ZZ + X and 
Z"f* + X, where dilepton pair masses are greater than 
30 GeV. This is the most restrictive cross section limit 
for ZZ production at the Tevatron. We set limits on 
anomalous neutral trilinear ZZZ and ZZ^j* gauge cou- 
plings. These represent the first bounds on these anoma- 
lous couplings from the Tevatron. Limits on ff , f§ , and 
f2 are more restrictive than those of the combined LEP 
experiments Q. 
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